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Abstract 

We extend previous techniques for calculations of transplanckian-energy string-string 
collisions to the high-energy scattering of massless closed strings from a stack of N Dp- 
branes in Minkowski spacetime. We show that an effective non-trivial metric emerges from 
the string scattering amplitudes by comparing them against the semiclassical dynamics of 
high-energy strings in the extremal p-brane background. By changing the energy, impact 
parameter and effective open string coupling A = gN, we are able to explore various 
interesting regimes and to reproduce classical expectations, including tidal-force excitations, 
even beyond the leading-eikonal approximation. 



1 Introduction 



One of the most remarkable properties of string theory is the generic presence, in its mass- 
less spectrum, of a spin-2 particle behaving as the graviton of general relativity. Invariance 
under general coordinate transformations emerges as part of the larger gauge symmetry 
of string theory and string interactions are automatically compatible with it. This is 
clearly reflected in the low-energy string effective action, which can be obtained from the 
scattering amplitudes of the massless excitations and which reduces to an extended - and 
typically supersymmetric - version of general relativity with additional short- distance 
corrections. As a result, in string theory spacetime is dynamical, even if the perturbative 
formulation of the theory requires the choice of a fixed background. Any highly massive 
or energetic state should then produce physical effects that can be ascribed to an effective 
spacetime metric. 

In a series of papers [1, 2, 3] 1 Amati, Ciafaloni and one of us (referred to hereafter as 
ACV) analysed the way in which classical gravity effects emerge from string computations 
performed in flat Minkowski space and developed methods to derive quantum and string 
corrections to general relativity. The process considered in ACV is the collision of two 
strings at trans-planckian energies in a flat Minkowski background. Given the strong 
gravitational fields generated in the process, one expects that the effective spacetime 
experienced by the two colliding quanta will be modified in a drastic way, possibly leading 
to gravitational collapse and black hole formation. 

In [1] it was shown that an eikonal resummation of the leading high-energy contri- 
butions to the S-matrix leads to results that are fully consistent with general relativity 
expectations whenever effects due to the finite string size (l s = y/a'h) can be neglected. 
In this particular process the effective geometry turns out to be, to leading order, the 
well-known Aichelburg-Sexl (AS) shock-wave metric [8] produced by an energetic point- 
like massless particle. This simple and intuitive interpretation unfortunately fails beyond 
the leading-eikonal (or small deflection angle) approximation [3]. In general the effective 
geometry produced by the collision is expected to be very complicated and to depend on 
the system as a whole rather than on each individual colliding string. An all-order ap- 
proach was proposed in [9] (see also [10]) but it turned out to be too difficult to implement 
this beyond the leading-eikonal approximation. 

Recently there has been some progress [11, 12, 13, 14, 15] in trying to compute these 
"classical corrections" to all orders through a dimensionally-reduced effective action, but 
the approximations made are not yet fully under control. Nonetheless these approxima- 
tions seem to lead to physically sensible results, at least in terms of identifying semiquan- 
titatively the critical parameters for gravitational collapse and of comparing them with 
those based on the formation of a trapped surface [16]. 

Several effects due to a non- vanishing string length l s ^ were also computed in [1, 
2, 17]: some of them do have a general relativity interpretation (e.g. as the effect of tidal 

1 See also [4] for a relevant string- loop calculation and [5], [6] for other derivations of the leading 
eikonal in the field theory limit. A different approach [7] is valid at fixed (and even large) angle within a 
high-energy window. The two methods agree in the intersection of their respective regions of validity. 
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forces exerted by an AS metric on extended objects [18]) while others (like the possible 
absorption of the elastic channel due to s-channel formation of heavy strings) do not. 
On the whole, a picture emerges whereby string-size effects prevent gravitational collapse 
when the Schwarzschild radius of the would-be back hole is smaller than the string length 
parameter l s while the approach to gravitational collapse is characterized, at the quantum 
level, by a rapid increase in multiplicity and by the corresponding softening of the final 
quanta [17, 11, 14]. The transition to the black- hole formation regime, which resembles a 
phase transition in general relativity, may turn out to be smoother in the quantum case. 

In this paper we shall apply the approach developed by ACV to the study of a different 
process, the scattering of a closed string from a stack of N parallel Dp-branes in Minkowski 
spacetime. The D-branes are massive solitons for which a microscopic string description is 
available [20]. This important property makes the string-brane system an ideal framework 
to understand the way in which string scattering amplitudes evaluated in flat space can 
provide information about the dynamics in an effective curved spacetime 2 . Indeed, from 
the point of view of perturbative string theory the presence of a collection of Dp-branes is 
entirely taken into account by the addition of an open string sector with suitable boundary 
conditions and does not require any modification of the background. On the other hand, 
from the point of view of the low-energy effective field theory the Dp-branes are a massive 
charged state and their presence will necessarily result in a curved spacetime. 

The backreaction of the D-brane system on spacetime is expected to be well-described 
by the extremal p-branes [21], which are BPS solutions of the supergravity equations of 
motion with a non-trivial metric, dilaton and Ramond-Ramond (p + l)-form potential. 
For p < 7 and in the string frame the extremal p-brane solution is given by 

ds 2 = 1 (r ]a pdx a dx 13 ) + y/H{r){5 ij dx i dx>) , (1.1) 
y/H(r) 

e«*> = g [H(r)}^ , C 01 ... p (x) = - 1 , (1.2) 

where the indices a, (3, . . . run along the Dp-brane world- volume, the indices . . . indi- 
cate the transverse directions and r 2 = 5ijX l x^ . Finally 

where g is the dimensionless string coupling constant, N the number of Dp-branes and 
Q n the volume of the n-dimensional unit sphere. This effective description should be 
reliable as long as the curvature is small in string units. Evidence that N parallel Dp- 
branes correspond to the curved spacetime given in Eqs. (1.1) and (1.2) was provided 
in [22, 23, 24] where it was shown that the large distance behaviour of the classical 
solutions can be recovered from string-brane scattering amplitudes at tree level. 



2 There is an analogue of this in quantum field theory: as shown long ago by Duff [19], a class of tree 
diagrams for the scattering of a test particle from a classical source reproduces the physical effects of the 
effective Schwarzschild metric generated by the source. The difference is that in string theory we have a 
microscopic quantum description of the source itself and of its couplings to the test particles. 
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We shall analyse the system of iV parallel Dp-branes using as a probe a high-energy 
massless closed string. In the high-energy limit the dominant interaction will be gravity 
and we will not be sensitive to the other background fields. We shall also work in the 
regime of small string coupling (5 < 1) and we shall therefore disregard all non-planar 
diagrams since their contribution is suppressed by powers of g. On the other hand, if 
we keep A = gN fixed (and thus let N — > 00) contributions with many open string 
loops (boundaries) are not suppressed. This is a welcome feature: indeed, in the regime 
where the energy E of the closed string is very large the scattering amplitude grows with 
energy in such a way that partial-wave unitarity breaks down at any finite order. These 
large contributions to the S-matrix from diagrams with a different number of boundaries 
are expected to combine together and yield a unitary eikonal form for the resummed 
amplitude. 

As in the case of the string-string scattering, our ultimate aim is to calculate the 
eikonal phase to all orders in the classical and string corrections. It will turn out that the 
problem of string scattering from a stack of D-branes is much easier to deal with than 
the original ACV-type collision problems. Indeed, the former problem is similar to that 
of scattering in an external field with the closed string acting as a probe and not back- 
reacting on the geometry. As a result, while the resummation of the closed string loops in 
the string-string scattering studied in ACV leads to an effective geometry that depends 
non-trivially upon the energies of the two colliding quanta, the metric that emerges from 
the resummation of open string loops considered in this paper (e.g the length R p defined 
in (1.3)) depends only on the D-brane system. 

Finally, we note that in our approach the existence of a non-trivial background is 
taken into account by the inclusion of surfaces with boundaries in the perturbative series, 
rather than by a direct modification of the couplings in the closed string sigma model. 
Since the string-brane scattering amplitude is evaluated using this microscopic definition 
of the system, in principle we can also analyse the background generated by N Dp-branes 
in regions where the effective description (1.1) and (1.2) is not reliable. 

The rest of the paper is organised as follows. In Section 2 we discuss the high-energy 
Regge limit of the two-point amplitude of a massless closed string state in the background 
of iV Dp-branes. We analyse in detail the disk and the annulus diagrams both in mo- 
mentum and in impact-parameter space and then, relying on the results of [1], we show 
that the resummed scattering amplitude has an operator eikonal form. We also compute 
explicitly the first subleading correction to the eikonal. In Section 3 we analyse the S- 
matrix by a semiclassical saddle point approximation at large impact parameters where 
one can expect that the description of the string-brane system in terms of a smooth ef- 
fective background metric should be recovered. In Section 4 we show that there is indeed 
full agreement with the results obtained by a semiclassical analysis of the propagation 
of a closed string in the extremal p-brane spacetime. In particular, we show that the 
string computation reproduces the classical deflection angle of a null geodesic to next- 
to-leading order as well as the inelastic amplitude for the string excitations due to the 
gravitational tidal forces. In Section 5 we discuss, with the help of a phase diagram, the 
different regimes that have been (or can be) explored in our approach and comment on 
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the relevance of this framework to address some aspects of black hole physics. Finally, 
in Appendix A and B we provide some more details about the string annulus calculation 
and the corresponding field theory diagrams discussed in Section 2. In Appendix C we 
do the same for the effective action (discussed in Section 4) that describes the excitation 
of a closed string moving in the curved spacetime sourced by the Dp-brane system. 

2 Scattering of a closed string on N Dp-branes 

Closed string two-point amplitudes in the background of N Dp-branes can be evaluated 
in perturbation theory summing over Riemann surfaces with an arbitrary number of 
boundaries and handles. In order to isolate the backreaction of the branes on the geometry 
of the spacetime one can consider the limit in which the open string coupling A = gN is 
kept fixed while g — > and N — > oo. The leading contribution then comes from surfaces 
with an arbitrary number of boundaries and no handles and the resummed amplitude is 
expected to allow an interpretation in terms of the propagation of a closed string in the 
curved spacetime of an extremal Dp-brane system, see Eq. (1.1). 

Our setup is reminiscent of the configuration considered in [25], where a highly ener- 
getic massless closed string was used to probe a collection of massive closed strings. The 
large number of constituents in the target plays in this case the same role as the number 
of D-branes in our case and can be used to suppress the contributions from the closed 
string loops with respect to the tree diagrams. Of course there are important conceptual 
and technical differences. For instance, a system of N parallel D-branes preserves 16 
supercharges; thus it cannot decay and the various constituents do not interact among 
themselves. This allows for a precise comparison between the results obtained from the 
string scattering amplitudes and those derived from the effective p-brane metric. 

In this paper we focus on the scattering of a massless closed string belonging to the 
NS-NS sector. The incoming and the outgoing strings are characterized, respectively, by 
their momenta pi and p 2 and their polarization tensors e\ and e 2 . The relevant vertex 
operators in the (—1, —1) and (0,0) pictures are 

7_ 1 (p,e) = fe p e-^ffe^ , 

V (p, e) = -^7^ (idX» + °^p^ (idX v + e tpX , (2.1) 

where k 2 = 2 & n 7 a' i g 2 is the gravitational coupling constant in ten dimensions and we 
follow the conventions of [26]. Translation invariance along the p + 1 directions parallel 
to the brane implies the conservation of energy and momentum along the D-brane world- 
volume 

(Pi+pa),, = . (2.2) 

The scattering amplitude is a function of two invariants of the residual Lorentz group 
SO(l,p) x 5*0(9 — p). Taking, without loss of generality, p\ and p 2 to have vanishing 
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components parallel to the brane, these can be written as 

s = E 2 = \ki\ 2 = \k 2 \ 2 , t = -(&! + A; 2 ) 2 = -4£ 2 sin 2 ^ , (2.3) 

where k = p± denotes the (9 — p)-dimensional vector transverse to the world-volume 
of the Dp-brane and G the angle between k\ and —k 2 . The two invariants correspond, 
respectively, to the energy of the incoming string and to the momentum transferred to 
the brane. For small deflection angles, k\ and k 2 are nearly (anti)parallel and their 
difference, {k\ — k 2 ), defines a privileged transverse direction, that of the large external 
momenta. The impact-parameter vector b, which lies in the direction of q = — (ki + k 2 ), 
is approximately orthogonal to such a direction and, together with it, defines the collision 
plane. In the following we shall often use bold characters to denote (8 — p)-dimensional 
vectors transverse both to the brane worldvolume and to the direction of (hi — k 2 ). 

Besides taking the large- N limit at fixed gN, we will also consider this process in the 
high-energy, small scattering angle limit defined by the conditions 

— > i , e < i , (2.4) 

h 

where we have re-introduced momentarily Planck's constant to emphasize that this limit 
allows for a semiclassical treatment. As already stressed in [1], this regime, even if 6 is kept 
fixed and therefore t becomes large, is dominated by soft dynamics, i.e. by the exchange of 
many gravitons whose individual transverse momenta are of order (and somewhat smaller 
than) Rp 1 <C \f—i and the corresponding relevant impact parameters b are somewhat 
larger than R p . This is the reason why we can justify the use of the Regge limit, order 
by order in the loop expansion, for the calculation of the scattering amplitude in b space. 
After resummation at fixed b, we can go back to t and verify that we can cover the 
kinematical region (2.4). 

Note that by fixing gN we are also fixing, through eq. (1.3), the ratio R p /l s . In the 
rest of this paper we will consider brane configurations with R p > l s , the stringy regime 
R p < l s will be briefly mentioned only in Section 5. When R p > l s the string-brane 
interactions are dominated by gravity effects and one expects a more transparent relation 
with the effective geometry in Eqs. (1.1), (1-2). Since the small- angle condition requires 
b ^> R p , we can see that our scattering amplitudes will be evaluated at impact parameters 
much larger than the string scale, b ^> l s . The conditions (2.4) also imply 

RnE bE J . . 

i"* 1, T = i ' (2 ' 5) 

where J is the (conserved) angular momentum. 

Finally, we will impose an upper limit on E such that the effects of closed string loops 
can be neglected. As we shall discuss more quantitatively at the end of this section, such 
a limit can be arbitrarily large if we go to sufficiently large iV (and thus sufficiently small 
g). When this condition is fulfilled, the two-point amplitude can be written as follows 

oo 

A(p 1 ,e 1 ;p 2 ,e 2 ) = i(2ir) p+1 8 p+1 (p u + m ) J^Mpu e l]P2 , e 2 ) , (2.6) 

h=l 
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where the label h counts the number of boundaries. Furthermore, in our limit the depen- 
dence of the scattering amplitude on the polarization tensors simplifies considerably, the 
leading term being proportional to Tr(eie^) 

A ft (pi,ei;p2,e 2 ) ~ Tr(e ie J) A h (s, t) , (2.7) 

where Ah contains only the part of the string amplitude that diverges at high energy. 
Setting 

oo 

A(s,t) = Y,A h (s,t) , (2.8) 
h=i 

we obtain the following expression for the Regge limit of the T-matrix, related to the 
S-matrix by S = 1 + iT, 

T( Pl ,e i; p 2 ,e 2 ) ~ (27rrV+V,|| Tr(e ie £) , (2-9) 

where the factor 2E takes into account the correct normalization of the asymptotic states. 
The first term in the series in Eq. (2.6) is the familiar disk amplitude, the second one is 
the annulus amplitude. We shall first analyse the Regge limit of these two amplitudes and 
then explain how to obtain the leading behaviour of the full perturbative series (2.6). As 
in [1] , the resummed amplitude takes a simple operator eikonal form in impact parameter 
space. The eikonal operator will allow us to study both the emergence of classical gravity 
effects from open string loops and the string corrections to the classical dynamics. 



2.1 Disk and annulus amplitudes 

The amplitude describing the elastic scattering of a massless NS-NS closed string state 
on the disk is well-known [22, 27, 28] 3 and reads 

A 1 {p 1 ,e 1 ;p2,e 2 ) = JC(pi, e x \ p 2 , e 2 ) - ^ _ — - ^ , (2.10) 

where the full kinematical factor /C(pi, ei; p 2 , e 2 ) can be found, for instance, in [22]. The 
amplitude exhibits poles both in the t-channel and in the s-channel. The former are due 
to the exchange of closed strings and appear when the two vertex operators approach 
each other, the latter are due to the exchange of open strings and appear when the two 
vertex operators approach the boundary of the disk. In the Regge limit the leading term 
in the kinematical factor is 

£(pi,ei;P2,C2) ~ Tr (eiej) (a's) 2 . (2.11) 

Even if in this paper we focus on states belonging to the NS-NS sector, we note that a 
similar result also holds for the elastic amplitude with two R-R massless fields. 

3 Disk amplitudes with both open and closed strings were first computed in Ref. [29]. 
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From (2.10) we obtain in the Regge limit 



Ms,t) = n -^0^T(-^t)e-^(a'sr^ . (2.12) 



-) V 4 



The previous formula shows that the amplitude is dominated by the exchange in the t- 
channel of the Regge trajectory of the graviton. The imaginary part takes into account 
inelastic processes where the closed string excites open string degrees of freedom attached 
to the brane worldvolume. As the energy increases the approximation of single-reggeon 
exchange eventually breaks down and the tree-level amplitude violates unitarity. As it is 
well-known both in field theory [30] and in string theory [1, 6], unitarity is recovered by 
taking into account multi-reggeon exchanges, described in our case by diagrams with a 
higher number of boundaries. 

Let us consider now the amplitude with two boundaries, the second term in the per- 
turbative series in Eq. (2.6). In general, perturbative string amplitudes are given in terms 
of correlation functions of the worldsheet theory integrated over the moduli space of Rie- 
mann surfaces with punctures. In the simple case of a surface with two boundaries and 
two punctures we have, besides the insertion points of the two vertex operators, Z\ and z 2 , 
a single purely imaginary modular parameter r op = ir 2 . Depending on whether the world- 
sheet time direction is chosen parallel or orthogonal to the two boundaries, this amplitude 
can be interpreted either as a one-loop diagram of open strings, the annulus diagram, or 
as a tree-level diagram of closed strings, the cylinder diagram. The two descriptions are 
connected by the modular transformation r op i— >■ t c \ — — 

As we will show in the following, at large impact parameters the dominant contribu- 
tions to the scattering amplitude come from the region of large A = Im(r c i). For this 
reason it is more convenient to represent the amplitude by means of two closed string 
vertices and two boundary states [31] and to display its explicit form in the closed string 
channel 

A 2 {p 1 ,€ 1 ;p2,e 2 )=Af J d 2 z 1 d 2 z 2 (B\V 1 (z 1 ,z 1 )V 2 (z 2 ,z 2 )D\B) , (2.13) 

where Af is a normalization factor and D is the closed string propagator. The final result 
is (see [32, 33] for a derivation of the amplitude and Appendix A for our conventions) 



A 2 (s,t) = 



o A 2 Jo Jo Jo Jo 



(2.14) 



In the previous formula the dependence on the external momenta is contained in the 
function 

I = e*^ , (2.15) 

which is the correlation function on the annulus of the exponential part of the vertex 
operators in Eq. (2.1). The functions V s and V t can be expressed in terms of the Jacobi 



7 



theta function 9\(z\t) (see Eq. (A.l) for the definition) and read 



2 , . fl 1 (iA(C + p)|<A)ei(<A(C-p)|<A) 



v s = + log ; 7'/ , (2.16) 



K = 8tt A Pl p 2 + log -7-^ ■ 77 77 77 , (2.17) 



0i(i\p + w 


i\)0i(i\p — uj 


iX) 


6>i(iAC + w 


iX)6 1 (iXC-uj 


iX) 



where we introduced the variables p = p\ — p 2 , C — P1+P2 and uj = uj\ — l>J2- The amplitude 
does not depend on a = ui 1 + w 2 as a consequence of its invariance under translations in 
the direction parallel to the boundaries of the cylinder. 



2.2 The Field Theory limit at high energies 



Before presenting the detailed derivation of the high-energy limit of the string amplitudes, 
which involves several technical points, it is useful to consider the disk and annulus di- 
agrams in the field theory limit. This limit, which corresponds to sending a 1 — > while 
keeping R p fixed, considerably facilitates the analysis and allows us to illustrate our main 
results in a simpler setting. The amplitudes discussed below can also be derived by evalu- 
ating directly the Feynman diagrams contributing to the scattering process, as described 
in Appendix B. 

The field theory limit of the disk amplitude is easily obtained from Eq. (2.12) and 
reads 



r(^) (-0 • 



(2.18) 



As expected, at high-energy the amplitude is dominated by the exchange of a single 
graviton in the t-channel. 

The field theory limit of the annulus diagram in (2.14) is more subtle since this ampli- 
tude is given in terms of an integral over the variables A, p, ( and 00. It turns out that the 
variable T = a'X plays the role of the Schwinger parameter of the field theory diagrams 
and thus has to be kept fixed in the limit. Since a' — > 0, this implies that the relevant 
region of integration is the region of large values of A. We can then rewrite Eq. (2.14) in 
the following way 



A 2 (s,t) 



3 2 R 1A - 2 P 



7T S 



r 2 (¥) 



2 1 ? 



dT 
1 2 



d( 



C ri 

dp doo e 
C Jo 



-a'sVs-^/Vt 



, (2.19) 



where now V s and Vt stand for the asymptotic form at large values of A of the functions 
in Eqs. (2.16) and (2.17), given in Eq. (A. 2). The first term in V s in (A. 2) shows that at 
high energy the integral over p is dominated by the region of integration around p ~ 
and therefore we can make the further approximation 



a'V s ~ 



-2nTp 2 



a'V t ~ -2ttTC(1 - C) • 



(2.20) 
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Substituting in Eq. (2.19), we see that the dependance on a' disappears from the annulus 
amplitude. We note that in this limit Eq. (2.19) agrees precisely with the field theory 
integrand for the first Feynman diagram in Figure 3 of Appendix B, evaluated around the 
saddle point y ~ in the Schwinger parametrization (see Eqs. (B.l) and (B.2)). 

In order to make the integral over p convergent we can continue it analytically to 
negative values of E 2 or equivalently we can rotate the integration contour in p from the 
positive real axis to the positive imaginary axis 

C dp e 2 ^ ~ — L- . (2.21) 
< V2TE 

Evaluating the two remaining integrals over T and ( we obtain 

4 , (M) _ ^il^i (rf r( £zi) ?mm , ( , 22) 



r 2 (^) 2 1 ? y/2E \\t\J v 2 ; r(6-p) 

where the upper index (3) indicates that this term diverges as E 3 . An alternative expres- 
sion for the previous amplitude, which will be useful in comparing it with the correspond- 
ing full string theory result, is the following 

.Met) 1 I ^ ( iE2R l-"^ Y 1 t = _, (2 23 , 

We see that A 2 can be written as the convolution in momentum space of two tree-level 
amplitudes A\. 

The individual perturbative amplitudes, like A\ and A 2 , are divergent for large values 
of E and therefore at sufficiently high energy they will violate unitarity. It turns out, 
however, that the leading contributions in energy coming from the disk, the annulus and 
the higher loop diagrams, corresponding to surfaces with more than two boundaries, are 
the terms of an exponential series and combine together to give just a phase in the S- 
matrix. This is most directly seen if we express the scattering amplitude as a function of 
the impact parameter b instead of t. The relation between the S-matrix and the amplitudes 
computed above is 4 

S=l + iT=l + i^ , (2.24) 
while momentum and impact-parameter space are related by a Fourier transform 



4 We omit here for simplicity of notation the terms containing the polarizations of the two external 
particles and the <5-function required by momentum conservation. The complete expression is given in 
Eq. (2.9). 
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Using the following equation 



/ 



d D k Ml ,, v v 2^ r(i/+f + 1) 1 



i e **(*T = -—o ^7, -77^ > ( 2 - 26 ) 



(27r)^ V ' Z/+f 7T^/ 2 r(l-I/) (fe2 

we can transform both Eqs. (2.18) and (2.22) to impact-parameter space, getting 5 




2b Q -p r(^=2) / ' ' 2 \ 2 6 6 -f r(^=2) 

In conclusion, we obtain 

S(E, b) = l + iT ± (E, b) + tT 2 (E, b) + --- = l + iT 1 (E,b)-± (T\(£, b)f + ... , (2.28) 

which is consistent with the fact that the unitarity violating disk amplitude exponentiates 
giving just a phase in the S- matrix 

S(E,b) =e m{E ' b) . (2.29) 

We can now proceed to compute the next-to-leading term in energy in the field theory 
limit. In this discussed in Appendix B, we need to take into account several 

Feynman diagrams. Their combination reproduces the a' — > limit of Eq. (A. 22) in 
Appendix A 



Z 

This result should be compared with the field theory integrand for diagram (a) in Figure 3 
of Appendix B, where we keep the subleading term in the expansion around the saddle 
point in the Schwinger parameter y. After performing the integrals over T and ( in Eq. 
(2.30) we obtain 



P(V)2' "/ V V « ; V - - , 

Using again Eq. (2.26), the previous expression becomes in impact parameter space 

iT a, (E b) = ; 4 2) _ _. ^er^ rgg) 

,T 2 (E,b) _ ,— - , ba _ 2p 4r(g _ p) . (2.32) 

The upper indices (2) and (1) denote that and diverge respectively as E 2 and E 
at high energy. In the next subsection we will generalize these results in order to include 
the string corrections. 



5 Here and many times in the following we encounter expressions which exhibit singularities at partic- 
ular values of p. The correct result is obtained by taking the limit from non-integer p, as in dimensional 
regularization, and in throwing away irrelevant infinities. As an example, at p = 6 eq. (2.27) produces 
logarithms of b besides an unobservable infinite Coulomb phase. 
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2.3 String corrections to the leading eikonal 

We now extend the results obtained in the previous subsection in two directions. On the 
one hand, we evaluate the high-energy limit of the full string amplitude (2.14) in order 
to derive the string corrections to the leading field theory result; on the other hand, we 
determine the form of the terms in A 2 that scale as E 2 , which are small in comparison to 
the leading behaviour E 3 , but still yield a divergent contributions to the T-matrix (2.9). 
In the previous Section we have already extracted those terms in the field theory limit, 
but their derivation presented in the Appendix includes also some (but possibly not all) 
string corrections. 

The high-energy limit of a string amplitude can be determined by evaluating the 
asymptotic behaviour of the multidimensional integral over the world-sheet moduli using 
saddle-point methods [34]. The leading terms are due to critical points in the interior or 
on the boundary of the integration domain, each critical point corresponding to a specific 
degeneration limit of a Riemann surface with punctures. We leave the technical steps 
necessary to evaluate the saddle-points of the annulus with two punctures in Eq. (2.14) 
to the Appendix A. Here we just give a sketch of the derivation and summarize the main 
results. 

In the limit (2.4) there are no critical points in the interior of the domain. The 
dominant contribution comes from the region of small values of p, since we are taking the 
high-energy limit for the external states, and from the region of large values of A, since 
we are considering impact parameters larger than the string scale. We shall then expand 
the integrand for large A and small p and perform exactly the integrals over the angular 
variables u>i and u 2 , which implement the level matching condition. Since we wish to 
focus on the part of the S-matrix (2.9) that is perturbatively divergent at high energies, 
we will consider only terms that are at least of order E 2 . It turns out that the expansion 
of the integrand around p = yields a leading contribution (indicated by A^) that scales 

as E 3 and a subleading contribution that scales as E 2 (indicated by A^). 

(3) 

The leading term A 2 in the annulus amplitude is precisely the one required for the 
exponentiation of the tree-level amplitude into an operator eikonal form, which generalizes 
the field theory result given in Eq. (2.23) and reconciles the perturbative string expansion 
with unitarity. More precisely, to all orders in a' the leading annulus contribution has the 
following form 6 (see Appendix A for details) 



6 Strictly speaking (2.33) is not well defined because there are non-physical divergences, due to the 
poles of the gamma function in V2, whenever k is such that a'(ti + £2 — t) = —2, —4, . . .. These diver- 
gences, however, arise from a region where the integrand is suppressed by factors of E. As shown in the 
Appendix A, the full amplitude is free of unphysical poles. 




(2.33) 



where 



v 2 {h,t 2 ,t) 



r[i + f(t! + t 2 -t)] 
r 2 [i + ^(t 1 + t 2 -t)} ' 



(2.34) 
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and 

q 2 = -t , h = -k 2 , t 2 = - (q - k) 2 . (2.35) 

The function V2(ti,t 2 ,t) is the vertex for the emission of two reggeized gravitons derived 
in [1] in the context of high-energy string-string collisions. The fact that precisely the same 
vertex appears in the string-brane scattering process studied in this paper is not surprising. 
As in the field theory description of the deflection of a particle by an external potential, one 
expects that the string will interact with the brane through multiple exchanges of states 
lying on the Regge trajectory of the graviton. The leading term of the annulus amplitude 
confirms this expectation showing that the amplitude factorizes into the product of the 
vertex for the emission of two reggeized gravitons and two disk amplitudes which encode 
the boundary conditions pertaining to the brane source. 

In the a't — > limit the string amplitude (2.33) reduces to the leading field theory re- 
sult (2.23) of Section 2.2 and a comparison between the two expressions helps to clarify the 
effect of the string corrections. We see that the graviton pole is replaced by the complete 
Regge trajectory of the graviton and that the convolution of the tree-level amplitudes 
now involves a non-trivial kernel V2(ti,t2,t). As it stands, it is not evident that (2.33) is 
the second term in the expansion of a simple eikonal form. In order to better understand 
the structure of the series let us consider the leading terms provided by the higher-order 
amplitudes. 

For a generic amplitude Ah with h boundaries, the term with the highest power of E 
should correspond to the exchange of h reggeons and scale as E h+1 . We can then write 

A h h+1 \s,t) 1 i h - 1 £=f f d^ki 



2E h\ {2E) h 



T-r f P k- 

II J J^=p Ms,ti)...A 1 (s,t h )V h (k 1 ,k 2 ,-M) , (2.36) 



where U = — k 2 for i — 1 . . . h and Yli=i ki = q generalize the conventions in (2.35), while 
the vertex for the emission of h reggeized gravitons Vh [1] generalizes the h = 2 expression 
in (2.34). To resum all the leading contributions in Eq. (2.36), it is important to realise 
that the vertices Vh have a simple representation in terms of vacuum expectation values 
of string vertex operators [1] 



h 2 - 



14(k 1; . . . , kfc) = <0| J] I g : : |0> , (2.37) 

i=i{ 

where the string fields X(a) are defined in (A. 17). By using this operator form for the 
vertex Vh, we can rewrite the integrand in (2.36) as a convolution in momentum space 
and then diagonalize the full string expression going to impact parameter space 



2E J (2tt) 8 -p 2E 



u 2tt 



1=1 n 
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Summing the contribution of all the /i-reggeon diagrams we finally obtain the leading 
high-energy behaviour of the series in Eq. (2.8) 



E 

h=i 



Ah j; ' j ~ (0|- \e 2i5 ^ - 1 |0) , (2.39) 
2£/ ! L J 



where 



2tt J (2tt) 8 -p 2£ ' 7 2tt 2£ 



(2.40) 

As anticipated, the eikonal form of the amplitude is the one expected in a high-energy 
process dominated by the exchange of reggeized gravitons between the string and the 
branes. The string eikonal differs from the field theory eikonal in that it involves the 
full string amplitude (2.12) and the impact parameter is shifted by the string position 
operator. It is this simple shift that takes into account effects due to the finite string size 
which can become relevant already at large impact parameters, as we shall explain in the 
next Section. 

We now turn to the subleading term A^ in the annulus amplitude (2.14). This 
subleading contribution has the same energy dependence as the tree-level amplitude (2.12) 
and its natural interpretation is that it induces a renormalization of the tree-level result 

A 1 (s,t) ^A (2 \s,t) = A 1 (s,t)+A 2 2 \s,t) + ... , (2.41) 

where the index (2) denotes that all the terms in the previous series diverge as E 2 at high 
energy. The correction to the Regge pole term remains non-trivial also in the field theory 
limit and it is crucial to reproduce the precise form of the Dp-brane solution in Eq. (1.1). 
For instance, in Section 4 we will compute the deflection angle of a null geodesic in the 
Dp-brane space-time and compare it with the string result in the limit of large impact 
parameters b. In this example one can see explicitly that the one-loop shift A 2 yields 
a non- vanishing contribution to the deflection angle which accounts for the classical sub- 
leading corrections in R p /b. 

The sub-leading contributions are computed in Appendix A. In particular in Eq. (A. 23) 
we give the complete expression for the term which is dominant in the a't — > limit. It 
is this term that, when multiplied with the prefactor in Eq. (2.14), reduces in the field 
theory limit to 

the expression already given in Eq. (2.31). 

Let us briefly comment on the subleading terms in the higher-order (many bound- 
ary) amplitudes. Besides the leading terms A^ +1 ^ which give the eikonal operator in 
Eq. (2.39), string diagrams with h boundaries should also contain contributions of type 
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A^\ that is with the same energy dependence as the tree- level amplitude (2.12). These 

(2) 

terms are similar to the term A 2 just discussed and should provide additional renor- 
malizations of the Regge pole term to be included in (2.41). Clearly there should also 
be several other terms scaling with a power of the energy intermediate between A^ +1 ^ 

(2) 

and A h . All these terms should combine to define a high-energy S-matrix which gener- 
alizes the leading eikonal operator (2.39), including both subleading classical corrections 
in powers of R p /b and string corrections in powers of a'/b 2 . Exponentiation of the former 
is necessary in order to reproduce the classical curved spacetime expectations discussed 
in Section 4. In order to determine the form of the latter one should perform a detailed 
analysis of the higher-order amplitudes which is beyond the scope of this investigation. 
Nonetheless, some indications about these string corrections to higher-order terms in the 
Rp/b expansion should be provided by generalizing our study of string propagation in the 
D-brane metric (1.1) presented in Section 4. 

To conclude this Section let us go back to the conditions under which we can neglect 
higher-genus topologies, i.e. closed string loops generating extra handles. Since closed 
strings live in the bulk any extra closed-string loop should involve the ten-dimensional 
Newton constant G w ~ k\ ~ g 2 (ct') 4 with no associated power of N. This is the reason 
we have already advocated to argue that, at sufficiently large N, closed-string loops can 
be neglected. In order to make the argument more quantitative let us notice that, in 
analogy with the ACV case, each extra handle can contribute to the eikonal phase at 
most one extra power of the asymptotic-energy variable, which is E in the present case. 
Hence, by dimensional arguments, each extra handle will be associated with a factor 
In general, from a surface with n boundaries and m handles we expect a contribution to 
the eikonal phase of the form 

Eb {R p \ n{7 - p) (G w E\ m Eb ( E\ m AT . n 2m fl s \ {7 ~ p)n+7m 



><■">"* ~T{T) l^J ~T[m) Wl-tf) ■ (2-43) 

The first correction to the disc amplitude (2.27) will thus be 

Eb ( R P \ 7 P f G\qE \ Eb E 3 r ( l x \ 14 p 



On the basis of (2.44) we see that the contribution of one extra handle is negligible if 

where / is a simple function of its arguments. This condition is obviously satisfied at 
fixed Rp/l s (i.e. fixed gN) and fixed Rp/b (i.e. fixed G) if iV is taken to be sufficiently 
large. It is easy to see that, under this condition, we also have S(n,m) ^ 1. Note finally 
that, given the above-mentioned relation between do and a', the condition (2.45) allows 
for E to be arbitrarily large even with respect to the 10-dimensional Planck mass. 
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3 Analysis of the amplitude 



The eikonal operator derived in the previous Section 

S(s, b) = e™^ , 2S(s, b) = -L J ^ : A, (a, b + X(<r)) : , (3.1) 

gives the leading behaviour, at high energy and in a series expansion in powers of R p /b, 
of the scattering amplitude of a massless closed string on a stack of N Dp-branes. This 
amplitude, which resums the dominant contributions of an infinite number of open string 
loops in Minkowski space, should be interpreted as describing the semiclassical propaga- 
tion of a closed string in a curved spacetime, the background generated by the D-branes. 

In this Section we shall use Eq. (3.1) to analyse the way in which this curved back- 
ground influences the motion of our probe string and to obtain some information on the 
background itself. In our approach we do not assume any knowledge of the effective metric 
generated by the D-branes, rather we see it emerging dynamically from the resummation 
of the open string loops. The fact that all-order perturbative string computations in flat 
space give rise to effects which can be interpreted as due to the motion of the string in an 
effective curved background was first discussed in the context of string-string collisions 
at transplanckian energies [1, 35]. As explained in the introduction, our results provide a 
simpler instance of the same phenomenon. 

Since the process we are considering is the scattering of an extended object in an 
asymptotically flat spacetime, there are two main effects we can study: the deflection of 
the trajectory of the projectile and the excitation of its internal degrees of freedom. As 
we will show, these two effects are neatly taken into account by the eikonal operator (3.1). 

In the next Section we will provide a direct derivation of the deflection angle and the 
excitation spectrum of a string in the extremal p-brane background of Eqs. (1.1), (1-2), 
finding precise agreement with the results of this Section. The agreement between the 
resummation of open string loops in Minkowski space and the quantization of the string 
in an external metric confirms that, for large R p , a collection of iV coincident Dp-branes 
is well approximated by the extremal p-brane solution of the supergravity equations of 
motion. It also shows that, at least in certain limits, our approach provides a quantitative 
tool to study the classical dynamics of a string in a curved spacetime. The advantage of the 
present approach, based on the microscopic and manifestly unitary D-brane description, 
is that it does not rely on the existence of an effective external metric and can be used to 
analyse other interesting dynamical regimes, as discussed in Section 5. 

In order to derive from (3.1) both the deflection angle and the excitation spectrum, 
we expand the eikonal phase in a power series in the string position operators X 1 



2 5( S ,b + X)~^ 



(3.2) 



assuming b ^> R p ^> l s ^\n(a's) and keeping only the first two non-trivial terms. The 
indices i, j label the 8 — p directions of the impact parameter space transverse to the 
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brane and to the collision axis and the symbol A denotes the average of a local operator 
A(a, t) on the worldsheet 



da : A(a,r = 0) : . (3.3) 



The two terms in (3.2) will give rise respectively to the leading contribution to the de- 
flection angle and to the leading contribution to the tidal excitation of the string modes. 
Let us analyse these two effects in turn. 

3.1 Deflection angle up to next-to-leading order 

The first term in (3.2) is simply the Fourier transform of the disk amplitude (2.12) 



where l s (s) is the effective string length, the size of a string of energy E = yfs 

l s (s) = l s Vlna's . (3.5) 

The previous formula shows that when b ^> R p and R p ^> l s (s) the eikonal phase is 
predominantly real, since the absorption effects due to the imaginary part in (3.4) becomes 
relevant only for b < l s (s). In the computation of the deflection angle p we can then 
approximate A\ with its real part. 

The functional relation between the deflection angle and the impact parameter can be 
derived evaluating the Fourier transform of Eq. (3.1) back to momentum space 

S(s,t) ~ J d 8 - p b e- ibq e ^ Al(s ^ . (3.6) 



The integral is dominated by the saddle point 

1 IdAi 



^ 2Eb db 

and, using the small-angle relation Q p ~ — we obtain 



b , (3.7) 



(2) 

In our analysis in Section 2, we showed that the annulus amplitude contains a term, A 2 , 
which has the same energy dependence as Ai but is of higher order in R p /b. To determine 
how these additional term enters in the eikonal operator in Eq. (3.1) would require a more 
detailed study of the subleading contributions coming from surfaces with more than two 
boundaries. It is however plausible, at least if one neglects the effect of string corrections, 
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that also this term will exponentiate and contribute to the series expansion of the eikonal 
phase in powers of R p /b. If we make this assumption, we can include in the real part of 
the eikonal phase in (3.8) the one-loop renormalization of the Regge pole obtaining 



^ r p=fc) r R^p 



Re[Al + A * )] ' ° v " r (^) tf-p J 4r(t£) r (5 - P ) b™-*> 



(3.9) 



Our string computation then leads to the following approximation for the deflection angle 
of a null geodesic in the background of iV Dp-branes 



EM ( R 
r(¥) 



7- P ir (15=2Ej 



2T(Q-p) \b 



Rr 



2(7-j>) 



+ 




3(7- P y 



(3.10) 



In the next Section we will show that this result is in perfect agreement with the deflection 
angle predicted by classical gravity in the extremal p-brane background, giving direct 
evidence in favour of our assumption that also the subleading term exponentiates. 



3.2 Tidal excitation at leading order 

We now turn to the excitations of the internal degrees of freedom of the string. They 
are taken into account by the second term in Eq. (3.2), which is of leading order both 
in Rp/b and in a' /b 2 . The higher-derivative terms in the Taylor expansion of the eikonal 
phase in Eq. (3.1) lead to higher string corrections weighted by the same leading power 
of R p /b but of higher order in a' /b 2 and are therefore suppressed when b 3> l s - It is the 
term with two derivatives that provides the leading contribution to the imaginary part of 
the eikonal phase, although it gives a negligible correction to its real part and therefore 
to the deflection angle. As we shall show it gives rise to absorption effects which, unlike 
those due to the imaginary part of the tree-level amplitude in Eq. (3.4), become important 
already at large distances b 3> l s (s) [1]. This imaginary part of the eikonal phase accounts 
for the excitation of the string under the influence of the long-range gravitational field of 
the brane [18]. 

In order to study this effect we derive the explicit form of the eikonal operator (3.1) in 
terms of the standard operators a l n , a l n) n e Z, corresponding to the left-moving and right- 
moving modes in the expansion (A. 17) of the fields X 1 . Since the eikonal phase depends 
only on the modulus of the impact parameter, the matrix of its second derivatives has 
the following simple structure 



1 d 2 A 1 (s,b) 



Q±(s,b) 



' bibj 



+ Q\\(s,b) ^, (3.11) 



with 7 — p coincident eigenvalues associated with the components of the string operator 
orthogonal to the impact parameter 

g ^ M) = v^^^' (3 - 12) 
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and one eigenvalue associated with the component parallel to it 

1 (fA^s.b) 

Q||M) = ^^^. (3.13) 

The average of the square of the string coordinates can be expressed in the following 
convenient form [1] 

(^) 2 = « J2 ~ [2^,n-l-n,n-7ln] , i=p + l,...,8, (3.14) 

1 

n=l 

where the operators T* are 

Ti, n = l^n^n , 27£ n = 1+1 (c/_ + , (3.15) 

and satisfy the following commutation relations 

\T-,m ^+,m] = 2 ^nm , Po,„, ^±,m] = ^ T±, n ^ij ^nm ■ (3.16) 

The previous relations imply that the operators T£ form a S77(l,l) algebra whenever 
they carry the same n and i labels, otherwise they commute. To quadratic order in the 
string oscillators the eikonal operator reads 

oo 8 

e 2iS(s, b} „ e ^W) Yl e i ^ Qfc ( 2T o*"- 1 - T +."- T -.») , (3.17) 

71=1 fc=p+l 

where stands for Q± or according to whether the direction k is perpendicular or 
parallel to the impact parameter. Using the identity 

e x(2T -T + -T_) = e -^T + e -21n(l-x)To ^ 

we can rewrite the eikonal operator in normal ordered form 

oo 8 ia 'Qk 

«*<-•>- n n fzis («•!») 

n=l fc=p+l n 

If e «-<«'«fc +-™e V n 7 l °' n ; e Q k -- n . 

n=l k=p+l 

It is now immediate to evaluate the elastic scattering amplitude, which is given by the 
vacuum expectation value of the eikonal operator 

(0|e 2 ^ (s ' b) |0) ~ e ^ Al{s ' b) e- ia '< {7 - p)Q±+Q ")r (l - ia'Q^) T 7 - p (1 - ia'Q ± ) , (3.20) 
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where we used the infinite-product representation of the gamma function and 7 is the 
Euler-Mascheroni constant. From the previous formula one can see that the string cor- 
rections to the real part of the eikonal phase are negligible at large impact parameters. 
Their most important effect in this limit is to induce a non-vanishing imaginary part 



which, after using the identity T(l + ix)T(l — ix) 



(0|e 



2i5(s,b) 



0} 



~ e 



. T x , can be written as follows 

sinn ttx ' 

o 1 

irrv'Ou 1 2 

(3.21) 



na'Qj_ 


7-p 
2 


Txa'Q\\ 


sinh 7ia'Q± 




sinh ira'Q\\ 



At high energy the previous formula becomes 



(0|e 2 ^ (s ' b) |0> 



-^7=Im^i(s,6) 



(2na')^\Q ± (s,b)\^\Q ll (s,b)\ 



e -|a'[(7-p)|Q x (s,6)| + |Q||( S ,6)| 

Using the explicit form of Ai(s, b) we obtain 



(3.22) 



Q±(s,b) = Q 1 (s,b) 



n^T (^) w-v 



2 v r (^=2) b 8 ~P ' 
The expression for the elastic amplitude simplifies 



Q^s,b) = -(7-p)Q 1 (s,b) . (3.23) 



2i5(s,b)\ 



| ) ^ e^^r (1 + io'(7 - p)Q 1 ) T 7 - p (1 - ia'Qi) , (3.24) 
(0|e 2 ^ (s ' b) |0) ~e-^ Im - 4l(s ' 6) (2W|g 1 (s,6)|)^ ^/y—p e -™'(7-p)\Q^)\ ? (3.25) 

and we can see that the absorption of the elastic channel due to string excitations becomes 
non negligible for b < bp where 



1 2 ) -Rl' p 



(3.26) 



At large distances this effect is more important than the inelastic absorption given by (3.4) 
and becomes relevant already at impact parameters large compared with both the string 
scale and the curvature of the brane background. In the next Section we will show that 
the tidal excitations of the string given by (3.25) agree with the results of a semiclassical 
computation in the Dp-brane spacetime. 



4 Comparison with curved-spacetime expectations 

In this Section we will compare the results we have obtained from string computations in 
Minkowski spacetime to what we would expect from the propagation of point or string-like 
objects in the non-trivial background (1.1), (1-2) generated by the D-branes. Note that, in 
the high-energy limit, our probes are only sensitive to the metric part of the background 
(since gravitational couplings are proportional to the energy), while the dilaton and the 
RR-form backgrounds would provide subleading corrections not considered in this paper. 
We will focus our attention on the comparison of two effects: 
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• Deflection angles up to next-to-leading order in the point-particle limit. 



• Tidal excitation of stringy probes at leading order. 

As we shall see, there is full agreement on both effects between the string-loop calculations 
and curved-spacetime expectations. 

4.1 Classical deflection up to next-to-leading order 

Let us start by a (quite standard) computation of the deflection suffered by a massless 
point-like probe in the metric (1.1) produced by our stack of Dp-branes. We shall use a 
slightly more general form of the metric allowing one to change conformal frame (e.g. from 
the string frame to the Einstein frame). Since the coordinates involved in the geodesies 
are only the time t and the spatial coordinates r and 8 of the plane in which the motion 
takes place, we can limit ourselves to the following part of the metric 

ds 2 = -a(r)dt 2 + f3(r)(dr 2 + r 2 d6 2 ) . (4.1) 

Since the metric does not depend on the time t and the angle 6 there are two conserved 
quantities, the energy E and the angular momentum J. From these two conservation laws 
and from the invariance of the action under arbitrary reparametrizations of the world-line 
coordinate u, one can find a differential equation relating 9 and the radial coordinate 

d9 b d9 b , in . 

~r = 7= & ~r= i — , 4.2 

dr ^yff^ dp vW--&V 

where b = J/E, p = R p /r, b = b/R p . Notice that the previous result depends only on 
the ratio a/ (3 and it is therefore invariant under an r-dependent rescaling of the whole 
metric. This means that both the string and the Einstein frame metric yield the same 
equation (4.2) for the classical trajectory. In the last step of (4.2) we used the actual form 
of a/ (3 given in Eq. (1.1). 

From (4.2) we see that the value of the angle 6 at the turning point r* is 

de f p * b 



6 ( r *) = / 1Z dr 



/ dp = , (4.3) 



/n ^Jl+p 7 -P-b 2 p 

where p* = R p /r* is the smallest root of the equation 1 + p 7 ~ p — b 2 p 2 = 0. Since the 
trajectory of a probe particle in the metric in Eq. (4.1) is symmetric around r*, the 
deflection angle Q p is given by 

Q p = 20(r,) - 7T => e p = 2 f P * dp b 7T . (4.4) 

Jo Jl+p 7 -P-b 2 p 2 
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The integral can be performed explicitly in terms of elementary functions for the cases 
p = 5, 6 yielding 



tan 



1 

2b 



7T 



— 7T 



(4.5) 



For the case p = 3 we get instead 

e 3 



where K is the complete elliptic integral of first kind and 



(4.6) 



h = -i + -[b- \b 2 -A 



(4.7) 



Similar expressions involving elliptic integrals can be given for p — 1 and p = 4 while 
for the cases p = and p = 2 we do not have an expression in terms of special func- 
tions. Nonetheless, the leading and next-to-leading terms in the large impact parameter 
expansion for arbitrary p can be computed and read 



35) 

r(^) 



Rr 



7-p ir (i5^ 



+ 



2T(6 -p) V b 



R, 



2(7-p) 



+ 



(4.8) 



in perfect agreement with the string calculations in Eq. (3.10). 

We should make a remark at this point about the order at which we expect to find 
agreement between the string calculation and the external metric one. In matching the 
two results we have identified the impact parameter b of the string calculation, defined 
by the Fourier transform (3.6), with the impact parameter of the geodesic calculation, 



defined as b = J/E. We expect such an identification to fail at order 
instance, sin© p starts to differ from Q p . 



when, for 



4.2 Tidal excitation of the closed string at leading order 

We shall now compare the results of the previous Section concerning the possible excitation 
of the probe closed strings with what one obtains by quantizing a closed string in the non- 
trivial metric (1.1). A similar exercise in the case of string-string collisions leads, to leading 
order, to agreement with expectations for quantizing a closed string in an Aichelburg-Sexl 
metric [36]. 

In the case at hand the end result of the string calculation is the eikonal-operator 
formula (3.1). Such a formula refers to the leading contribution in R p /b but is supposed 
to hold at all orders in a'/b 2 . On the other hand, the curved-spacetime calculation we 
shall present below is limited to small string fluctuations around the point-particle null 
geodesic while, in principle, it can be extended to higher orders in Rp/b. Our comparison 
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will be made in the overlap of the domains of validity of the two calculations, namely at 
leading order both in a'/b 2 and in R p /b. In spite of this the perfect agreement between 
the two calculations appears to be almost miraculous, given the very different techniques 
being used, and represents in our opinion a highly non trivial check of the validity of our 
approach. 

In order to set up the curved space calculation we shall follow [37] and rewrite the full 
metric as 

ds 2 = a{r) i^-dt 2 + J2(dx a ) 2 ^J + (5{r) (dr 2 + r 2 {d9 2 + sin 2 9dQ 2 _ p )) , (4.9) 

where in our case f3(r) = l/a(r) = y/H{r). This reduces to (4.1) on the plane of the 
null geodesic considered in Sect. 4.1 and we are now interested in describing the metric 
around such a geodesic. This is done by first going to a system of adapted coordinates 
u, v, z, x a , y l in which the geodesic corresponds to constant v, z, x a , y 1 and u = u(r) plays 
the role of the affine parameter along the geodesic 

dv = -dt + bd9 + Cdr , dz = d{9 + 9(u)) , (4.10) 



du = ±—, C(r) = J^-j--. 4.11 

C V air) r2 



Here 9{u) is the angular coordinate 9, evaluated along the null geodesic as in (4.2) and 
expressed in terms of u via Eq. (4.11). The physical meaning of the 7 coordinates x a , y % is 
that they represent fluctuations orthogonal to the null geodesic and, respectively, parallel 
to the directions of the brane world- volume (x a ) or along the (7 — p) directions (y l ) which 
are orthogonal both to the brane and to the plane of the geodesic. The z coordinate is 
orthogonal to the brane, but lies in the plane of the geodesic. In our conventions the 
point u — corresponds to the turning point r* and the choice of sign in the equation 
for u depends on the point of the geodesic we are considering: we choose the minus sign 
in the approaching region and so parametrize the part of the geodesic from infinity to r* 
with the interval — oo < u < 0; for the remaining part from r* to infinity we choose the 
plus sign and so it corresponds to the interval < u < oo. In these adapted coordinates 
the metric takes the form 

ds 2 = 2dudv - adv 2 + Ibadvdz + r 2 aC 2 dz 2 + adx a dx a + f3r 2 sin 2 (z - 9)dQ 2 _ p .(AA2) 

At this point we can take the Penrose limit of the above metric (corresponding to the high- 
energy limit for the probe) and focus on a small neighborhood around the null geodesic 
by expanding to the quadratic order the dependence on all coordinates transverse to the 
light-cone and keeping only the linear terms in v. This clearly eliminates the dv 2 and 
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dvdz terms, then by the following change of coordinates 
-.0 
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(4.13) 
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we can bring the metric in the pp-wave form [37] 

V 7-p 

ds 2 = idudv + y d %l + Y dy2 + d y° + g( u > v\ y°) du2 , 
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The bosonic part of the string sigma model then reads 



(4.14) 



(4.15) 



where from now on we will drop the hats on the coordinates. In the previous equation So 
is the free Minkowski string action, r] a p the flat worldsheet metric and Q is as in Eq. (4.14) 
but now considered as a function of the string coordinates (denoted by capital letters). 
In these coordinates string quantization is quite easy if we choose (within the class of 
orthonormal gauges) the light-cone gauge 7 

U(a, t) = ap u T -> a'Er . (4.16) 

This choice drastically simplifies the non-trivial part of the sigma model action to give 

e r 2n da r +oc 

S-S = - — dug(u,X a (a,u/a'E),Y t (a,u/a'E),Y°(a,u/a'E)) 

2 Jo 27T 



E r 2 * da 
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du Q x {u) Y 0) + Q y {u) Y Y?{°, 0) + Q (u)Y 2 (a, 0) 
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F r 2w rlrr ( P 7 ~ P \ 

\ a=l i=l / 



(4.17) 



where in the second step we have used the high-energy limit. In this way the integrals 
over u decouple from the string coordinates and just provide c- number coefficients Q to 

7 The factor 2 usually present in the r.h.s of Eq. (4.16) is absent here because we take < a < 2n. 
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the quadratic action of the fluctuations. Because of the change of sign in du/dr at u = 
the integrals under consideration are twice the same integrals between and oo. 

At first sight all the fluctuations of the closed strings, both in the Neumann and in 
the Dirichlet directions, appear to be excited. It turns out, however, that the former are 
not excited to the leading order in R p /b to which we are working. This is because we can 
write 



■f 

Jo 



-du 



/ (d u In Va) 2 + d 2 u (in ^a) du . (4. 18) 

Jo L - 1 



, \2(7-p) 

In the second integral the first term is proportional to ( -y 1 and so it yields a 

contribution of higher order in R p /b, while the second term, being the total derivative of 
a function that vanishes on the integration boundaries, gives zero (see Appendix C for 
more details). Hence, to this order, c x = 0. 

This is not true for the other two fluctuations. The coefficient of the Y-i fluctuations 
is quite simple to evaluate by writing, using again the trick that we have just used, 
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d u In y//3r sin 6^ + d 2 u (in ^/(3r sin 9} 



du . (4.19) 

( R V~ p 

As in the previous case, the first term does not contribute at the order ( -y 1 we are 

interested in, but, as shown in Appendix C, this time the u = boundary provides a 
non- vanishing contribution for the second term and we obtain c y = —Q p /b. 

Inserting the leading value of P from (4.8) we see that the coefficient of Y 2 agrees 
with the one obtained in Section 3 from the second derivative of the phase shift in the 
directions orthogonal to the brane and to b, i.e. 

7-p -n(8-p\ 
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Computation of the coefficient of the Yq fluctuations is a little more involved since the 
analogue of the square of the first derivative in (4.18) is not subleading: one has to perform 
carefully the second derivative, expand it to leading order in R p /b, and integrate explicitly 
each term over u. This lengthy but straightforward exercise (reported in Appendix C) 
gives for the coefficient of the F fluctuations 



co = ^(7 - p) ) ^ =► f co = Q„ (*, b) , (4.21) 

again in perfect agreement with the string calculation. 

5 Conclusions 

The eikonal operator derived in this paper allows us to describe different regimes of the 
scattering of a closed string off a stack of N Dp-branes. These regimes are characterized 
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by the relative magnitude of the scales involved in the dynamics, namely the energy of 
the colliding string E, the impact parameter b, the curvature radius R p of the brane 
background and the string length l s . In analogy with ACV, it is useful to draw, as in 
Figure 1, the different regimes in an (R p , b) plane marking on both axis the effective string 
length l s (s) = ZsVTn a's. Notice that, in our case, the parameter R p , which sets the scale 
of the effective geometry, is independent of E and we can therefore describe the various 
regimes at a fixed large energy a's 3> 1. We have also limited the diagram to the region 
R p , b > l s (s) since, in this paper, we have not discussed the new phenomena that occur 
for b or R p smaller than l s (s). 




Figure 1: Various qualitative regimes in the scattering of a closed string of fixed high 
energy off a stack of Dp-branes as a function, in a log-log plot, of R p and b, both taken 
to be bigger than the effective string scale l s (s). The different regions are discussed in the 
text. 

There are various distinct regimes that we shall now describe by going from the top- 
left to the bottom-right part of Figure 1. For very large b (above the blue line) the eikonal 
phase is small and we are in a perturbative - rather than in a classical - regime. This 
region corresponds to infinitesimal deflection angles (as b becomes larger and larger). The 
straight black lines at 45 degrees represent curves of constant deflection angle: the solid 
one is meant to represent a critical ratio of R p /b below which, classically, a test point 
particle is captured, while the dashed one represents a typical, small deflection angle 
9 C 1. The red line represents instead the impact parameter bp, defined in Eq. (3.26), 
below which tidal excitation becomes relevant. 

In the region bordered by the blue, the red and the black-dashed line the leading 
eikonal approximation is reliable and elastic unitarity holds to a very good approximation. 
Below the red and above the black-dashed line we have the tidal excitations computed in 
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Section 2.3, while in the complementary region (above the red but below the black-dashed 
line) we need to use an improved eikonal including higher-order classical corrections whose 
first term was computed in Section 2.2. Tidal effects, however, should be small in this 
region. Clearly, as we go further down in the diagram, we encounter a region in which both 
kinds of correction (string-size and classical) come together. When the impact parameter 
is close to R p , we expect large corrections to the eikonal operator discussed in this paper 
and, for this reason, the solid line, signaling the onset of the string-tidal effects, ends 
before meeting the b ~ R p line. 

In this paper we have focused on the upper part of the diagram in Figure 1, where the 
impact parameter is large compared to R p and gravity effects dominate the interaction be- 
tween the Dp-branes and the string probe. It is in this region that the comparison between 
the string dynamics in the extremal p-brane background and the dynamics resulting from 
the string scattering amplitudes is most transparent. Having tested our approach in this 
regime, there are two other interesting but more difficult regions to consider. The first 
is the string region R p < Z s (s), where string corrections to the geometry are important 
and one expects that the dynamics will be very different from the one predicted by the 
effective background in Eqs. (1.1), (1-2). The second is the region b ~ R p > l s (s), where 
the dynamics should be dominated by strong gravity effects. As we lower the impact pa- 
rameter to study these two new regions, we should also be able to make contact with the 
analysis of high-energy amplitudes at fixed angle as discussed in [7] and, in the context 
of D-branes, in [38, 39]. 

The quantitative analysis of a string-brane scattering process in the string region R p < 
l s (s) and in the strong gravity region b ~ R p requires some control on both the classical 
and string corrections to the leading eikonal operator. An example of the possible effects 
of these corrections is provided by the string excitations along the Dp-brane world-volume. 
Although these excitations are absent in (2.40), the analysis of Section 4 suggests that 
they should become relevant at smaller impact parameters. An effective way of studying 
these kinds of corrections is to calculate amplitudes similar to those discussed in this paper 
but involving also massive string states [40]. From a more general perspective, it would 
be very interesting to derive and interpret the corrections to the eikonal operator which 
are of a higher order in ^ and by computing explicitly the high-energy behaviour 
of amplitudes with three or more boundaries. This kind of analysis might be possible by 
using Regge-Gribov techniques [3]. 

Another interesting generalization of our setup is to change the nature of the mas- 
sive target and possibly also the nature of the asymptotic space by including compact 
directions. For instance, bound states of different D-branes, such as the D1/D5 system, 
have been studied intensively in recent years and provide a tractable system that has an 
exponential degeneracy of states (see for instance [41] and references therein). The ge- 
ometries corresponding to these states are known and have some interesting features such 
as the breaking of rotational symmetry in the transverse space. It would be interesting to 
generalize the string computations of [42, 43] and to see how the properties of these mi- 
crostate solutions affect the high-energy scattering considered in this paper. Furthermore, 
by adding a momentum charge to the D1/D5 bound state, one obtains the configurations 
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studied in [44] which are related to a black hole with a macroscopic horizon. 

In all these examples the most interesting region from the point of view of gravity 
and black- hole physics is that of small impact parameters (b < R p ), where, as previously 
noted, the eikonal phase receives large classical corrections. However, as was the case for 
the simpler regime analysed in this paper, the mixed open/closed string amplitudes which 
describe the collision of a string with a D-brane configuration might be more tractable at 
high energies. They could then provide an ideal framework for discussing in a quantitative 
way the dynamics of matter falling beyond a horizon and the resulting small perturbations 
of a black hole above its ground state. 

Acknowledgements 

We would like to thank Daniele Amati, Costas Bachas, William Black, Marcello 
Ciafaloni, Emilian Dudas, Stefano Giusto, Vishnu Jejjala, Lello Marotta, Cristina Monni, 
Yaron Oz and Gabriele Travaglini for fruitful discussions. This work is partially supported 
by INFN and by STFC under the Rolling Grant ST/G000565/1. GD would like to thank 
the members of the Laboratoire de Physique Theorique de l'Ecole Normale Superieure in 
Paris, where part of this work was done, for the warm hospitality extended to him. 

A The saddle point in the annulus moduli space 

The string amplitude (2.14) contains an integral over the moduli space of a cylinder 
with two punctures in its interior representing the external closed string states. Our 
parametrization of this surface is summarized in figure 2. 




Figure 2: The string world-sheet is represented by the annulus between the outer circle of 
radius 1 and the inner circle of radius e~ wX . The two external states are represented by 
the two punctures located at Z\ and z 2 . 

We are interested in a particular kinematics where both the impact parameter and 
the energy of the incident states are much larger than the scale fixed by the string length: 
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ya'E 3> 1 and b 3> l s (s). In this regime the integral over the world-sheet moduli is 
dominated by the contribution of a critical region at small values of p = p\ — p<i and large 
values of A. Indeed, the large b limit implies that the momentum exchanged between the 
string state and the D-branes is small and therefore the corner of the string integration 
region that is more relevant is A — > oo. Moreover, we are taking the high-energy limit for 
the external states and, as we will see, the integrand is exponentially suppressed unless 
p ~ 0. 

The Jacobi theta function 

oo 

e^r) = -2e mT/4 sin(7rz/) ]~] [(1 - e 2mmr )(l - 2 cos(27rz/)e 2 ™ + e 4 ™)] , (A.l) 

m=l 

which appears in the definitions (2.16) and (2.17) has a simple expansion in our critical 
region. We can then derive the following expressions for the large-A behaviour of V s and 

V t 

V s ~ -2vrAp 2 - 4 (sin 2 + sinh 2 TrAp) (e' 2 ^ + e - 27rA(1 ~ c) ) , (A.2) 
V t ~ -2ttA [C(l - C) + P 2 ] + In [4 sin 2 iru + 4 sinh 2 n\p] . 

In the expansion for V t we kept only the contributions that do not have any exponential 
factors, while for V s we also included the exponential terms, since they are sizable as soon 
as A > men's. The first term in each expression corresponds to the leading field theory 
contribution considered in Section 2.2, while the other terms yield the sub-leading field 
theory contributions and the string corrections. 

Let us now focus on the integral in the square parenthesis of (2.14). The integrand 
does not depend on a = u>i + u 2 and is periodic under the shift u> — > u> + 1, so we can 
perform the integral over a and obtain 

/»oo 1*1 /*1 

/ = / d\ / d( / dp dco A^ X . (A.3) 
Jo Jo J~c Jo 

Since the integrand is an even function of p, we can restrict the integral over p to the 
interval from to (. The leading contribution in the large E limit scales as E 3 and 
comes from the region where p is set to zero, except for the first term of V s , which has 
a different structure and has to be kept exact. The subleading contribution is related to 
the expansion of X at the quadratic order around p = and we will see that it scales as 
E, i.e. with two power less of energy with respect to the leading contribution. Higher 
orders in the p expansion are further suppressed in the large E limit. So we can further 
approximate the building blocks (A.2) of the integrand and write X ~ X;X S with 

J,(p) ~ e w ^ 2+27rA « 1 -^^ , (A.4) 
Up) ~ (4sin 2 7 r W + 4 7 r 2 AV)-^e 4a ' s (- 2 ^^V)(e--c +e --( 1 -o ) ^ 
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where we wrote explicitly only the dependence on p. We then separate the two contribu- 
tions mentioned above by writing I — I± + I 2 

h = / d\ d( dp dul l (p)l s (p = 0) , (A.6) 
Jo Jo Jp~o Jo 

h= / d\ d( dp dwl^p) X s (p)-l s (p = Q) . (A.7) 

Jo Jo </p~0 Jo 1 J 

Let us first describe how to evaluate the leading term Ji, following closely [1]. We can 
regularize the integral over p by performing a Wick rotation on the energy of the external 
states E — > iE e . Then we perform the saddle point integration over p and, after Wick 
rotating back to Minkowskian energy, we obtain 

h = r d\ C d( f 1 dcuX^e^-^Up = 0) . (A.8) 

V2a's Jo Jo Jo 

In order to deal with the integral over uj, we expand the exponential in X s (p = 0) as 
follows 



OO OO 



n\m\ 

n=0 m=0 



and perform the integral by setting x = ttoj and using 

nn/2 

B{a,b) = 2 / dx sin 2 "" 1 x cos 2 ^ 1 x . (A. 10) 

Jo 

The integral over A can then be evaluated using the integral representation of the gamma 
function and the result reads 

Ji = i ^ y f 1 ^ 1 (4^r-r(^)i?(l + n + m-f,|) 

V^s^Jo ^! 7r ^ + i[2nC + 2m(l-C)-C(l-C)f]^ ' 



In this equation ( plays the role of a Feynman parameter for a diagram with two propa- 
gators and an integrated momentum over 8 — p dimensions 



I 



(2ra')^r(^) f d s -Pk 1 



o [2n( + 2m(l - C) - C(l - Of ] ^ 7 (27r)§ " t 2n + ^ \- 2m + ^ k " ' 

(A.12) 

where q is a (8— p)-dimensional vector whose norm is q 2 = —t. In terms of this momentum 
integral, the contribution l\ becomes 



1 ~'v / 2^(27ra') f ~ 4 n t^ 7 ( 27r ) 8 ~ P [2n + f k 2 ] [2m + f (k - q) 2 ] ' 

(A.13) 
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We can now transform the sums over n and m into two integrals over the complex plane 
by using 

1 f(m)s m 



E 

m=0 



m\ m + t 



dz m f{Zm)Sm e -inZr, 



c 2ni z m + t 



(A.14) 



and a similar expression for the sum over n. The contour C runs anti-clockwise around the 
real positive axis in the complex plane (clearly this identity holds if / does not have any 
pole on the real positive axis). Deforming the integration over the two complex variables 



we pick the contributions of the poles on the negative real axis 
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(A.15) 
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^l^4(-k 2 + k.q), 2 



(4e- l V S )^ k2+k ^| , 



where we included only the first pole in the Euler beta function in (A. 13), since the 
remaining poles (in z n + z m + 1/2 — a't/4 = —k with k > 1) yield contributions that 
are finite or small in the high-energy limit. Each line in this equation contains some 
non-physical poles when k takes particular values. For instance, the first term of the 
integral is divergent when 2 — a'(t + k 2 ) = 0, but this singularity is cancelled by a similar 
pole in the second line of (A.15). Similarly the first pole of the beta function in the last 
line is cancelled by the other pole of the second term (while the remaining poles of the 
beta function would be cancelled by the terms we neglected in (A.15)). By including the 
prefactor in (2.14) and using the duplication formula T(x)T(x + |) = \fii 2 1 ~ 2x T(2x), we 
can finally rewrite the last term of (A.15) as in (2.33). 

The Euler beta function in the last line of (A.15) can be written in terms of the 
expectation value of two tachyon-like string vertices [1] 



2tt 



i=l 



-1 : : |0> = —B - + ¥ k lk2 , - 



(A.16) 



where X(z) has the standard mode expansion, but without the zero-modes describing the 
centre of mass position and momentum 



E 



a n c inc _|_ 



n 



n 



(A.17) 



which yields the OPE X(z)X(w) ~ —a' In |1 — w/z\. At this point one can simply follow 
the same steps as in Section 2.3 and use (A.16) with k x = k and k 2 = q — k in order to 
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show that which follows from the last line of (A. 15), can be written as the convolution 
in Eq. (2.38). 

Let us briefly comment on the leading string corrections at small a't (or at large impact 
parameters b), which are contained in the last line of (A. 15). They are obtained by first 
expanding for small a't the r.h.s. of (A. 16), or equivalently the vertex V 2 in (2.34), 

V 2 (k,q-k)~l + y ^yk.(q-k)y + ... , (A.18) 

and then by using this approximate form in (2.33). One can then see explicitly that also 

the contribution of these first string corrections takes the form of a convolution involv- 
es) 

ing the second derivatives of A\ as in (3.11). After rewriting A 2 in impact-parameter 
space, we can compare this result with the expansion of (3.21) for small a'Q± and small 
a!Q\\ finding agreement at the quadratic order. The result in (3.21) shows how these 
string corrections are properly incorporated in the eikonal operator once diagrams with 
an arbitrary number of boundaries are resummed. 

We now analyse the subleading contribution I 2 . For simplicity we shall focus on the 
term that is non- vanishing in the field theory limit, that is the first term in the series 
expansion of the exponential in the definition of X s in (A. 5). The other terms can be 
evaluated in a similar way. 

With the change of variable sin 2 ixuj — > x, one can see that the integral over oj yields 
a hypergeometric function 

/ du (sin 2 ttu! + y 2 ) c = V 2c 2F 1 (-c, 1/2; 1; -y' 2 ) 
Jo 

1 r(| + c) yi+acW-i-c) 

; + y v 2 ; + . . . , (A.19) 



ttIYI + c) r(-c 



where in the second line we wrote the first two terms in the small-?/ expansion. The first 
term in (A.19) is cancelled by the contribution of X s (0) in the definition (A. 7) and thus 
we obtain 

h ~ SwS^S^ ^P) 1 '^^ ' (s+l > 2+27rA C(1 -°- • (A.20) 



o p~0 



The expansion in the second line of (A.19) is appropriate for the regime of very high 
energies at large, but fixed, values of the impact parameter; in this way pA, which is of 
the order of b/(a'E), is small. If we keep increasing b and keep E fixed, we reach a point 
where y in (A.19) is of order one and so it is more appropriate to expand the result (A.19) 
for small c = — a't/ 4 at arbitrary values of y. By using for instance Appendix B.l of [45], 
we obtain 

y 2c 2 Fi(-c, 1/2; 1; -y' 2 ) = 1 - 2cln2 + 2c7i\p + 0(c 2 ) , (A.21) 

where in this case we used the full expression for y = sinh(7rAp), since \p is not necessarily 
small. The p-independent terms coincide with the small c expansion of the first term in the 
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second line of (A. 19) and so they have been already included in I±. The term proportional 
to 7i\p in (A. 21) agrees with the small c expansion of the second term in (A. 19) (remember 
that in that equation y is approximated with ir\p). This shows that the first order in 
t of Eq. (A. 20) captures correctly the contribution of I 2 also at very large values of the 
impact parameter b > a'E. 

At high energies, the integral over p in Eq. (A. 20) is dominated by the saddle point 
p ~ and after a Wick rotation E — > iE e we obtain the following asymptotic behaviour 
for I 2 in the s 3> t limit 



(A.22) 

The higher powers of p in the expansion of the hypergeometric function in (A. 19), which 
we have neglected, would yield contributions that are suppressed by further powers of 
1/E. This is a consequence of the relation between the expansion in powers of p and the 
asymptotic dependence on the energy mentioned after Eq. (A. 3). The integrals over ( 
and A are straightforward and yield 

/2 ^ (as) 2^f-r(i + ft) ( ., /} v - r O~~ * 

(A.23) 

From this equation, including the prefactor in (2.14) and sending a' — > 0, one can derive 
Eq. (2.31). 



B Dilaton scattering in field theory. 

In this Appendix we provide some details about the field theory diagrams that contribute 
to the dilaton scattering from a stack of Dp-branes. As in the rest of this paper, we are 
interested in the high-energy limit of the full scattering amplitude and so we can focus on 
a subclass of diagrams. At the first order in (R p /b) 7 ~ p there is just a single diagram with 
a graviton exchange between the probe dilaton and the Dp-branes. It is straightforward 
to check that the Feynman rules summarized in Table 1 yield the high-energy result of 
Eq. (2.18). 

We can use the same Feynman rules to compute the diagrams in Figure 3. They scale 
as (R p /b) 2( - 7 ~ p " > and should therefore contribute both to the exponentiation of the leading 
term, see Eq. (2.22), and to the subleading term, see Eq. (2.31). We also need to consider 
a second class of diagrams, depicted in Figure 4, which are of the same order in R p /b 
and involve the exchange of gravitons, dilatons and RR-fields. While diagram 3(a) gives 
both leading (~ E 2 ) and sub-leading (~ E) contributions to the S-matrix, diagram 3(b) 
and the three diagrams in Fig. 4 give only sub-leading contributions. It is interesting to 
note that diagram 3(b) cancels, in the high energy limit and in the De Donder gauge, the 
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Graviton propagator 



- tn 1 




Dilaton propagator 
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2ik [k 1(jl k 2l/ ) - Ihpk^rj^] 




i(2n) 2 ^kx^^pT]^) - \{k 1{p k 2u) r] pu + k 1{p k 2a) r] pu )] + ... 









Table 1: Here we collect the Feynman rules necessary for the derivation of the diagrams 
in Figure 3. For the graviton propagator we use the De Bonder gauge. As usual, the 
round parentheses imply a symmetrization of the enclosed indices: for instance we have 
k\( [i k.2 V ) = (ki p k 2u +ki u k 2p ) /2 . The symbol r] pu means that the Lorentz indices are restricted 
to thep+1 directions along the brane world-volume. The dots in the quartic vertex stands 
for terms that are proportional to k\ ■ k 2 which cannot give contributions growing with E 
in the diagrams we are interested in. 



diagrams in Fig. 4. As a result, both the leading and the sub-leading contributions to the 
S-matrix come from diagram 3(a) alone. 



Figure 3: Diagram (a) contributes both to the leading (~ E 2 ) and to the subleading (~ E) 
terms in the S-matrix. Diagram (b) yields contributions scaling at most as E. 

Let us start our analysis from diagram (3) (a); by using the Feynman rules we have 



where pi, p 2 are the momenta of the two external states and k is the momentum flowing 
in the dilaton propagator. Notice that the combinations (pi — k) and (p 2 + k) are non- 
trivial only in the (9 — p) transverse directions, while k has also a time component and 
so k 2 = —E 2 + k\. By performing the same Wick rotation [E = iE e ) used in the 
string computation, we get standard Euclidean propagators and can introduce Schwinger 
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Figure 4: The contributions of these diagrams to the S-matrix scale at most as E. The 
third diagram involves the exchange of a RR-field represented by the dashed lines. At high 
energies they combine and cancel exactly the leading term of diagram 3(b). 



parameters for evaluating (B.l). By focusing on the square parenthesis we have 
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(B.2) 



In the high-energy limit, this integral is dominated by the saddle point y ~ 0. This is 
the same pattern seen in the string computation of Appendix A and indeed y plays the 
same role as p. In order to derive the leading term, one can neglect y in the second 
exponential and treat the integral over y in the Gaussian approximation. By inserting the 
result in (B.l) and using (2.6) we obtain the full result (2.22) for the dilaton amplitude 
at the next-to-leading order in the energy. 

The subleading contribution is obtained by expanding to linear order in y the second 
exponential in (B.2). Again by focusing on the saddle point y ~ all integrals can be 
performed and we obtain the contribution from diagram 3(a) to A 



(2) 



r(g) 
r(6-p) 



(B.3) 



(2) 

which agrees with the amplitude A 2 of (2.31). 
Diagram 3(b) yields 



3(b) = UKT v f{2^Y +l ^ +1 \p l+ p 2 ) 



d g - p q E 2 {7-p) 
(27if-Pq\q-{p 1 +p 2 )Y 

9-p 



(B.4) 



where q is the momentum of one of the two graviton propagators and the dots stand for 
subleading terms in E. The leading energy result in (B.4) cancels the analogue contri- 
bution from the Feynman diagrams of Figure 4. The direct evaluation of these diagrams 
is somewhat involved, in particular for the contribution containing the 3-graviton vertex. 
However, as an intermediate check, one can combine the three diagrams, before sewing 
the common dilaton-dilaton-graviton vertex, and see that, after a Fourier transformation, 
this reproduces the Dp-brane metric in the Einstein frame at the second order in the 
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(R p /r) 7 ~ p expansion. After including the vertex with the external states, we obtain, at 
leading order in E, a contribution opposite to (B.4). This shows that, in the gauge we 
chose, all contributions growing with the energy are captured by diagram 3(a). It would 
be very interesting to check whether a similar pattern holds also at higher orders in R p /b 
and to see whether the only relevant contributions at high energy are the "half-ladder" 
diagrams which are a natural generalization of 3(a). 



C Computation of cq, c x and c. 



Let us start by showing that c x in (4.17) is vanishing at order 



7-p 



As mentioned in 



the main text, the first term in the square parenthesis in (4.18) is clearly of higher order, 
since hiy/a ~ 0(R p /r) 7 ~' p . The second term is a total derivative and its contribution 



a. In 



d r In y/a 



0, 



(C.l) 



where we first used (4.11) to change the derivative with respect to u into a derivative with 
respect to r and then used C(r*) = 0. 

A similar computation holds also for c y . The combination r sin# = b + O ((R/r) 7 ~ p ) 

and so again the first term on the r.h.s. of (4.19) can be neglected at the order (j^j 
At leading order we have 



7-P 



C„ ~ 



2<9 u ln[ v //3rsin#(w)] 



d r \n[\ffir sin0(r)] 



(C.2) 



The function C(r) vanishes at the turning point r*, so the only possibility for obtaining 
a non-zero result is to act with the derivative on the sin 9 term. From (4.2) we see that 
this yields a factor of 1/C and the lower extremum of (C.2) yields the non- vanishing 
contribution 

6cos(^) 



bcos9 



r 2 /3 sin 9 



Q, 



rl/3(r*) sin 



Bp 



(C.3) 



where in the first step we used the relation between ^(r*) and G p in Eq. (4.4) and then 



we kept only the terms of order 



7-p 



. Using (4.8), we obtain the result (4.20). 



Let us finally consider the coefficient cq. From Eq. (4.14) we have 



c 



b 2 a 



\J j3r 2 — b 2 a 



du , 



which we can rewrite as 



c 



dr 



a/ /3r 2 — b 2 a 



d r d u \/ (3r 2 — b 2 a . 



(C.4) 



(C.5) 
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We then compute 



d u y/Pr 2 -b 2 a = ^d r ^(3r 2 -b 2 a = a' 1 ' 2 (\ + 



r JL _ h —°L 
2/3 2r~j3 



(C.6) 



where we have used Eq. (4.11) and a prime denotes differentiation with respect to r. 

(r \ 7 ~ p 

Since we are interested only in the lowest order terms that behave as ( J , we can 
expand a and j3 at large r as follows 



a ~ 1 + C a 



R r 



7-P 



R, 



7-P 



(C.7) 



We keep the coefficients C a and Cp arbitrary at this stage, but we will see that the result 
depends only on the difference Cp — C a , which is equal to one both in the string and in the 
Einstein frame. Substituting Eqs. (C.7) in (C.6) and taking the derivative with respect 
to r we get 



8 r 



R, 



7-p 



P-7 fR. 
2r 



7-P 



-C a + (p-7)Cp--(p-7)C c 



-C a + Cp(p-7)--C a (p-9) 



(C.8) 



We can now insert this expression in Eq. (C.5) and make the approximations, valid at the 
lowest order that we are considering, r* ~ b in the lower extremum of integration and a, 
(3 ~ 1 in the square root in the denominator. Using the integral 



dr 



•1-27 0f fcl-27 r( 7 



b 2 



r( 7 ) 



(C.9) 



we obtain 



C = f (7-P) 



R r 



7-P 



(C.IO) 



Since Cp — C a = 1, the previous result reproduces Eq. (4.21). 
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